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A field theory is presented for particles which interact via Coulomb and hard-core potentials. We 
apply the method to the one-component plasma (OCP) with hard cores, consisting of identical par- 
ticles of fixed charge and diameter in a neutralizing background, and the symmetric two-component 
plasma (TCP) with hard cores, consisting of equal numbers of positively and negatively charged 
particles of identical size. We obtain exactly the first few coefficients of a systematic low-density 
r*-) • expansion of the free energy for both models. The OCP coefficients go over to the classical Abe 

result in the high-charge (or small hard-core diameter) limit. The TCP, on the other hand, exhibits 
diverging coefficients in the high-charge limit, which is due to the formation of strongly bound 
ion-pairs. 
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PACS. 61.25.Hq, 83.70.Hq, 61.41. +e 



I. INTRODUCTION 



The equation of state of a one-component plasma (OCP), i.e. of a system consisting of point-like charged particles, 
is of central importance for plasma theory Q. It is commonly assumed that the OCP follows from the two-component 
plasma (TCP), consisting of a globally neutral mixture of positive and negative point particles, if one charge species 
is of lesser charge and forms a neutralizing background. Much interest in the plasma community has been focused on 
the limit of high plasma density (or large plasma parameter), which defines the limit when the electrostatic interaction 
at average separation between two particles is much larger than thermal energy H. This situation is, for example, 
realized in white dwarfs and the interior of large planets (such as Jupiter). Our interest in plasmas comes mainly 
from the close analogy with colloidal and electrolyte solutions, which also consist of charged particles. Very often 
r»o ' in colloidal systems, the concentration of charged species is rather low, such that a low-density expansion becomes 
meaningful. Also, in addition to Coulombic interactions, the excluded-volume interaction becomes important in 
colloidal and electrolyte systems, and therefore has to be included p|. The OCP without hard core interactions has 
been studied by a variety of methods II). The first few terms in a low-density expansion of the free energy have 
been obtained by an explicit resummation of the Mayer expansion H. The OCP with hard cores has in the past 
been studied extensively with numerical and integral-equation methods [pi . For the TCP, the leading term in the free 
energy has been calculated by Debye and Huckel in their classical paper [pj. The formal foundation for the rigorous 
thermodynamic treatment of a TCP including a hard-core interaction has been given by Mayer J7[] and Edwards || , 
but to our knowledge, no systematic calculation of the low density expansion of the free energy has been performed. 

In this article, we present a novel field-theoretic method which can be used to systematically calculate free energies 
and other quantities of interest for systems of charged hard-core particles in powers of particle densities. We apply this 
method both to the OCP and TCP in the presence of finite hard-core radii. In the limit of vanishing particle charges, 
we recover the virial expansion of a hard core gas for both the OCP and TCP. In the limit of vanishing hard-core 
radius, we recover the classical Abe result for the OCP Q|. The expansion coefficients of the TCP diverges in the 
limit of vanishing hard-core radius, due to the proliferation of bound particle pairs. Much interest has been focused 
on the critical behavior of the TCP with hard core interactions(the so-called restricted primitive model, RPM) |j] [l3[ . 
Our results demonstrate that the critical point of the RPM is not accessible by a low-density expansion. 

In the next two sections we present our methods and calculations for the TCP and OCP with hard cores, respectively. 
Our main results consist of the density expansion coefficients of the free energy, which are obtained in closed form for 
the full range of coupling constants. The final section is devoted to a brief discussion of our results. 

II. TWO-COMPONENT PLASMA 

We consider a gas of N + positive charges and AT_ negative charges, interacting with Coulomb interactions v(r) = 
Ib/t and an additional short-ranged potential w(r) which we assume to be a hard-core interaction of range a. The 



Bjcrrum length £b = q 2 /^neksT measures the distance at which the interaction between two charged particles equals 
the thermal energy. The canonical partition function of the system reads 
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is the charge density. The infinite particle self-energies v(0) and w(0) are subtracted. We perform a double 
Stratanovitch-Hubbard transform and get 
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where Z v and Z w denote the square root of the determinants of the Coulomb and hard-core operators, Z v ~ \/det v 
and Z w ~ vdetw, respectively. Assuming that the fugacitics of the positive and negative ions are the same, the 
grand-canonical partition function is defined by 
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and reads 
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where we used the definition 



h(r) = 0-^(^+^(0)72^ 
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We next introduce the Debye-Hiickel propagator, whose inverse is given by the ordinary screened Coulomb potential 
i)DH(r) = £bc~ kt /t with k 2 = 8tt£b^- The partition function can be written as 
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where 



Q(r) = h(r)e v{0)/2 cos[0(r)] - 1 - w(0)/2 + ^ 2 (r)/2 (7) 

and the average brackets denote averages with respect to the DH propagator and the hard-core propagator. The 
determinant of the Debye-Hiickel interaction is denoted as Zdh ~ Vdct wdh and will be explicitly calculated below. 
Defining the effective interaction Q in this fashion allows us to perform a systematic low-density expansion, since all 



higher cumulant expectation values of Q are finite even if the fugacity inside the Debye-Hiickel propagator is set to 
zero. The logarithm of the partition function is given by a cumulant expansion in powers of Q, 



\nZx 

V 



= 2A + Au(0) + V- 1 In 



Zdh 

~z7 



2AZi + 2X 2 Z 2 



where 
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All expectation values appearing in these expressions can be calculated explicitly. We note that, neglecting the hard- 
core interaction, the cumulant terms Z n correspond to giant cluster graphs considered by Abe m . Since the integrals 
with respect to the fluctuating <f> and ip fields decouple, we do not have to consider mixed expectation values. Up to 
second order in Q we only need the expectation values 
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(0 2 (n) cos[0(r a )]) = e-^(°)/ 2 (vdh(0) - t£ H (n " 'a)) • 
Introducing the difference between the Debye-Hiickel and the bare Coulomb self energy, which is finite, 

Av Q = v(0) - ^dh(O) = £bV8tt£b\ 
the functions Z\ and Z 2 can be explicitly written as 

Zx = c Av °/ 2 - 1 - Aw /2, 
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The DH contribution to the free energy in eq.(|h follows as 



Au(0) + V _1 ln 



Zdh 



dq ln f 1+ 8^A )+A 



(2tt 



dq 4tt£b (8tt£ b X) 



3/2 



(2tt) 3 g 2 



12tt 



(20) 



(21) 



Next we have to systematically carry out the back-Legendre transform to the canonical ensemble. To do this, we 
expand the grand-canonical potential in eq.(H) in fractional powers of the rescaled fugacity, A = a 3 A, 
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From the definition of the grand-canonical partition function, eq.fQ), it follows that the the rescaled concentration of 
ion pairs, c = a 3 (N + + iV_)/2V, is given by 
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From eqs.(22) and (E3|) we obtain the fugacity as a function of the concentration 
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Inserting this into eq.(E2), the rescaled free energy density follows as 
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with the expansion coefficients given by 
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To carry out the expansion of the grand-canonical free energy shown in eq. (g2|) , we also have to expand the functions 
Z n in fractional powers of the fugacity, 



Introducing the energy scale e 
obtain from eq.(plj) 



Z 2 = 4 0) + A 1 / 2 z«+0(A). 
IbI a with a being the hard-core radius (or diameter of the spherical particles), we 
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From eq.(20) we obtain 
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where Shi is the hyperbolic sine-integral function defined by 
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The coefficient of the A 3 / 2 term in eq.(p2|) follows from Eq.(pT|) as 
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the other coefficients are 
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In the limit e — ► 0, which corresponds to the pure hard-core case, we obtain hi = +87r/3, the well-known second- 
virial result. In the strong-coupling limit, e — ► oo, we obtain 62 = — 47re c /e. It is seen that the coefficient 62 
diverges exponentially as e increases. We attribute this divergence to the formation of strongly bound ion pairs: 
the Boltzmann weight of a pair of oppositely charged ions which stick together is e e , which dominates the statistics 
for large values of e. The next-leading coefficient 6 5 / 2 is dominated by its logarithmic contribution and is given by 
b 5 / 2 = — 2(27r) 3 / 2 e 9 / 2 ln(e) both in the limit of small and large e. In Fig.l we plot the free energy coefficients d 2 and 
d 5 / 2 as a function of e. For moderately large values of e, it is seen that the coefficients increase exponentially and 
thus render the low-density expansion useless for large couplings and large densities. It follows that the critical point 
of the TCP, which is supposed to occur at coupling strengths of the order of e « 15 and a rescaled ion pair density of 
c ss 0.03 |H]], is inaccessible by such a low-density expansion. 

As we will demonstrate in the following section, the OCP does not contain ion pairing and thus the expansion of 
the free energy in powers of the density is well-behaved even in the limit e — > 00. 



III. ONE-COMPONENT PLASMA 



We now consider a gas of N identical charged particles, interacting with Coulomb and hard core potentials. We 
assume the existence of a uniform oppositely charged background, which exactly neutralizes the charges. At first, we 
will neglect this background charge, and identify and omit the (infinite) energy contribution of this background in 
our final free energy expression. The canonical partition function of the system reads: 
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where p(r) = 2»=i ^( r — r «) denotes the particle density. The infinite particle self-energies are substracted. As we 
did for the TCP, we perform a double Stratanovitch-Hubbard transform and get 
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The grand canonical partition function is defined by 
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and reads 
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with the same definition for h(r) as in eq.(^). As for the TCP, the Debye-Hiickel propagator is defined by its inverse 
«DH(r) = ^Be~ Kr /r, but the screening length in this case reads k 2 = 4tt£b^- The partition function can be written as 
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where 



Q(r) = ^(r)e" (0)/2 -^ (r) - 1 - w(0)/2 + 2 (r)/2 
The logarithm of the partition function again follows by cumulant expansion, 
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with Z\ and Z 2 defined as in eqs.(|9j) and (10). The difference between the Debye-Hiickel and the bare Coulomb self 
energy reads 



Av = «(0) - udh(0) = 2i B y/nl B X, 



the functions Z\ and i?2 are given by 
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The DH contribution to the free energy follows as 



Aw(o) + v- 1 In 



Zdu 

Z v 



(4tt£ b A) 3 / 2 



12tt 



(43) 



(44) 



Next we have to systematically carry out the back-Legendre transform to the canonical ensemble. To do this, we 
expand the grand-canonical potential in fractional powers of the rescaled fugacity, 



^ = A- 63/2 A 3 / 2 -6 2 A 2 ... 

The rescaled concentration of ion pairs, c = a 3 N/V, is given by 

Xa 3 d\nZ x /V 

~ c = dx ' 

The fugacity as a function of the concentration is given by 
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The free energy density follows as 
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The prefactors are given by 
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The leading term of Z 2 is 
2?ra 3 



v + ^-(e 3 [7H-lne + r(e)-r(3V / 4^A)-ll/6] - e~ e [2 - e + e 2 ]) + 0(A) 



(51) 



and the leading term of Z\ 



Z 1 = ^+0(A 2 ). (52) 



The infinite constant v = J drv(r) corresponds to the interaction of the ions with the oppositely charged background 
and the background self energy and will be neglected in the following. The quadratic coefficient of the free energy 
reads 

d 2 = -|(e 3 [7 + lne + r(e)-r(3V4^)-ll/6J - e" e [2-e + e 2 ]) . (53) 

In Fig. 2a we plot this coefficient as a function of e in the limit of small coupling strengths and setting c = 1 inside 
the logarithm. As one can see, in the zero-coupling limit, e — > 0, this coefficient becomes 

2tt , N 

d 2 = y, (54) 

the standard second- virial result. We observe a non-monotonic behavior of the coefficient d 2 , and for e > 3.5 it 
becomes in fact negative. 

For large values of e, which corresponds to strong Coulomb interaction or relatively small hard-core diameter, the 
free energy density is conveniently rescaled by the Bjerrum volume £ B and reads 

e%f = cine + d 3/2 c 3/2 + d 2 c 2 + ■■■ (55) 

where we introduced the rescaled concentration c = £ B c (which is related to the so-called plasma parameter T, which 
is commonly used to characterize a OCP, by T 3 = 47rc/3). The coefficient appearing in this free energy expansion are 

2 

d 3 /2 = -gV^) ( 56 ) 

cfc = ~ (7 + he + r(e) - r(3 ^4^7?) - 11/6 - e'V' [2 - e + e 2 ] ) . (57) 

In Fig. 2b we plot d 2 as a function of e for large values of e, again setting the concentration appearing inside the 
logarithm equal to unity. In the strong-coupling limit, e — > 00, the coefficient approaches the value 
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which agrees exactly with the result obtained by Abe, by performing an infinite resummation of terms of the Mayer 
expansion |4J| . His calculation was done without hard-cores, and thus constitutes a special case of our more general 
result for arbitrary hard-core diameter. In Fig. 2b we denote the Abe limit, d 2 — — ^ [27— 11/6 + ln(3-\/47r)] — 
— 1.76476 by a broken line. As one can see, the convergence towards the asymptotic Abe result is quite slow, and even 
for strongly charged spheres with e = 10 the actual coefficient d 2 has reached only one half of its asymptotic value. 

IV. DISCUSSION 

In this paper we introduce a novel field-theoretic formulation for plasmas and charged electrolyte or colloidal 
solutions. We are able to include, in addition to the Coulomb interactions, hard-core or excluded volume interactions. 
We apply our method to the one-component plasma (OCP) and the symmetric two-component plasma (TCP) and 
obtain, as a main result, the first few terms in a systematic low-density expansion of the free energy. The major 
advantage over previous calculational methods is that each order in the density expansion, which corresponds to an 
infinite Mayer sum, is reexpressed as a single diagram involving a composite, non-linear, but local operator. This 
entails that our results are non-perturbative with respect to the coupling strength e = £e/a, and are thus valid 
both for vanishing electric charge, e — > 0, in which case we recover the virial expansion for a hard-core gas, and for 
vanishing hard-core diameter, e — > 00. In the strong-coupling limit, e — ► 00, our results for the OCP asymptotically 
approach the previous results by Abe, but even for values of the coupling parameter as large as e ~ 10 the differences 
to the asymptotic result are quite large. The TCP, on the other hand, does not have a well-defined strong-coupling 



limit, and we find the expansion coefficients of the free energy to increase progressively as the coupling parameter 
increases. This fact is connected with the formation of strongly associated ion pairs. One of the motivations of the 
present work was to elucidate the critical behavior of the TCP, which shows a critical point at a coupling strength of 
roughly e ~ 15. It is clear that unless resummation techniques of the low-density expansion are used, the low-density 
expansion cannot be used to predict the location of this critical point. 



[1] 
[2] 
[3] 
[4] 
[5] 

[6] 

[7] 

[«] 

[9] 

[10] 

[11] 
[12] 

[13] 



M. Baus and J.-P. Hansen, Phys. Rep. 59, 1 (1980); P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987). 

C. Deutsch, Y. Furutani, and M.M. Combert, Phys. Rep. 69, 85 (1981); S. Ichimaru, Rev. Mod. Phys. 54, 1017 (1982). 

I. Borukhov, D. Andelman, and H. Orland, Phys. Rev. Lett. 79, 435 (1997). 

R. Abe, Prog. Theor. Phys. 22, 213 (1959); E.G.D. Cohen and T.J. Murphy, Phys. Fluids 12, 1404 (1969). 

E. Waisman and J.L. Lebowitz, J. Chem. Phys. 56, 3086 (1972); F. Lado, Mol. Phys. 31, 1117 (1976); J.-P. Hansen and 

J.J. Weis, Mol. Phys. 33, 1379 (1977); D. MacGowan, J. Phys. C 16, 59 (1983). 

P.W. Debye and E. Hiickel, Z. Phys. 24, 185 (1923). 

J.E. Mayer, J. Chem. Phys. 18, 1426 (1950). 

S.F. Edwards, Philos. Mag. 4, 1171 (1959). 

G.R. Stell, K.C. Wu, and B. Larsen, Phys. Rev. Lett. 37, 1369 (1976). 

M.E. Fisher and Y. Levin, Phys. Rev. Lett. 71, 3826 (1993). 

J.M.H. Levelt Sengers and J.A. Given, Mol. Phys. 80, 899 (1993). 

N.V. Brilliantov, C. Bagnuls, and C. Bervillier, Phys. Lett. A 245, 274 (1998). 

R.R. Netz and H. Orland, Europhys. Lett., in press. 



FIG. 1. Plot of the coefficients di and ^5/2 in a systematic and exact low-density expansion of the free energy for the 
symmetric two-component plasma (TCP). The coefficients are shown as a function of the coupling strength e = Ib/cl, where 
£b — q 2 /(^neksT) is the distance at which two charged particles interact with thermal energy, and a is the hard-core diameter. 
It is seen that the coefficients diverge exponentially as e increases. 

FIG. 2. Plot of the quadratic free energy coefficient in the small-coupling limit, a), and in the strong coupling limit, b). For 
e — > 0, one recovers the hard-core virial coefficient, d,2 — 27r/3, and for e — > 00 the Abe result is obtained, which is denoted by 
a broken line. 
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Fig. 1c, Netz and Orland 



d 



o 



-40 



5/2 



-80 



-120 



0.2 




Fig. 1d, Netz and Orland 
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Fig. 2a, Netz and Orland 
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Fig. 2b, Netz and Orland 
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